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PROCEDE D'ESTIMATION DE DISTANCE POUR UN MOBILE SOUMIS A 
DES CONTRA1NTES DYNAMIQUES DE PARCOURS 

L'invention concerne la navigation de terrain pour un mobile 
soumis.a des contraintes de parcours variant au cours du temps, tel qu'un 
aeronef limite en taux de montee, la limite pouvant etre negative, et evoluant 
au-dessus d'une zone de terrain presentant des reliefs ou des obstacles 
menagants proches ou superieurs a son altitude de vol. 

Divers systemes ont ete developpes pour prevenir I'equipage d'un 
aeronef d'un risque de collision avec le sol. Certains, tels que les systemes 
TAWS (acronyme de I'expression anglo-saxonne : "Terrain Awareness and 
Warning System"), font une prevision .de trajectoire a court terme pour 
I'aeronef a partir des informations de vol (position, cap, orientation et 
amplitude du vecteur vitesse) fournies par les equipements du bord, la 
placent en situation par rapport a une carte de la region survolee extraite 
d'une base de donnees d'elevation du terrain accessible du bord et emettent; 
des alarmes a destination de I'equipage de I'a&ronef a chaque fois que la, 
trajectoire previsible a court terme entre en collision avec le sol. Ces; 
systemes TAWS agrementent leurs alarmes de recommandations.!. 
rudimentaires du genre "Terrain Ahead, Pull up". Certains d'entre eux, : 
donnent egalement des informations sur le niveau de risque de collision que • . 
font encourir les reliefs et les obstacles envirohnant I'aeronef sous la forme ft 
d'une carte presentant les reliefs ou les obstacles du terrain survole en V 
strates de couleurs differentes. Cependant, cette carte de risques de collision 
avec I'environnement tient uniquement compte des altitudes du relief 
relativement a la position du mobile et ne prend pas en compte de I'existence 
ou non d'une trajectoire realiste permettant de rejoindre les zones affichees. 

Pour satisfaire ce besoin de connaTtre les points du terrain survole 
restant accessibles apres une manoeuvre d'§vitement d'un relief ou d'un 
obstacle au sol, la carte de risque de collision avec I'environnement doit 
afficher uniquement les zones pour lesquelles il existe un chemin possible 
depuis la position courante du mobile. La realisation d'un tel affichage passe 
par I'association d'une metrique a une carte du relief tiree d'une base de 
donnees d'elevation du terrain. 

Une methode connue pour associer une metrique a une carte du 
relief tiree d'une base de donnees d'elevation du terrain a maillage regulier 
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de la surface terrestre ou d'une partie de celle-ci, consiste a considerer la 
carte presentant le relief a partir de valeurs d'altitude figurant, avec les 
coordonnees geographiques : latitude et longitude des points de mesure 
dans les elements de la base de donnees d'elevation du terrain comme une 
image dont les pixels sont les valeurs d'altitude des points de la base de 
donnees d'elevation du terrain illustrees dans la carte avec pour 
coordonnees abscisse et ordonnee au sein de I'image, les coordonnees 
geographiques latitude et longitude de ces points figurant dans les elements 
de la base de donnees d'elevation du terrain et a faire appel a une 
transforms de distance operant par propagation pour estimer des distances 
au sein de cette image. 

Les transformees de distance operant par propagation egalement 
connues sous I'appellation "transformees distance de chanfrein" ("chamfer 
distance transform" ou "chamfer euclidean distance transform" en anglo- 
saxon) deduisent la distance d'un pixel dit pixel but par rapport a un autre 
pixel dit pixel source, des distances precedemment estimees pour les pixels 
de son voisinage, par un balayage des pixels de I'image. Le balayage permet 
d estimer la distance d'un nouveau pixel but par rapport au pixel source par 
recherche du trajet de longueur minimale allant du nouveau pixel but au pixel 
source en passant par un pixel intermediate de son voisinage dont la 
d.stance a deja ete estimee, la distance du nouveau pixel but a un pixel 
intermediate de son voisinage dont la distance a deja ete estimee etant 
donnee par application d'un masque de voisinage communement appele 
masque de chanfrein. 

Une transforms de distance de ce genre a ete proposee en 1986 
par Gunilla Borgefors pour estimer des distances entre objets dans une 
image numerique, dans un article intitule :" Distance Transformation in Digital 
Images." et paru dans la revue "Computer Vision, Graphics and Image 
Processing", Vol. 34 pp. 344-378. Un des interets de ces transformees de 
distance par propagation est de reduire la complexity des calculs d'une 
estimat.on de distance en autorisant I'emploi de nombres entiers. 

Pour selectionner le trajet de longueur minimale donnant 
leshmation distance, une transforms de distance par propagation doit tester 
tous les trajets possibles. Cette obligation se traduit par une contrainte de 
regularity imposee a I'ordre de balayage des pixels d'une image G 
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Borgefors propose, pour satisfaire cette contrainte de regularite, de balayer 
les pixels d'une image deux fois consecutivement, dans deux ordres inverses 
I'un de I'autre, qui sont soit I'ordre lexicographique, I'image etant analysee de 
gauche a drorte ligne par ligne et de haut en bas, et I'ordre lexicographique 
inverse, soit I'ordre lexicographique transpose, I'image ayant subi une 
rotation de 90°, et I'ordre lexicographique transpose inverse. II propose 
egalement d'adopter un masque de chanfrein de dimensions 3x3 avec deux 
valeurs (3, 4) de distances de voisinage ou de dimensions 5x5 avec trois 
valeurs (5, 7, 1 1) de distances de voisinage. 

Les transformees de distance operant par propagation sont deja 
employees dans le domaine de la navigation terrain pour robot. Dans ce 
cadre, il est connu d'utiliser la transformee de distance de G. Borgefors avec 
une contrainte statique consistant a attribuer, de maniere autoritaire, une 
distance infinie a un point en analyse lorsqu'il apparaTt qu'il fait partie de. 
reliefs ou d'obstacles a contourner repertories dans une memoire de zones ,. 
interdites de franchissement, cela de maniere a eliminer, de I'ensemble des 
trajets testes lors d'une estimation de distance, ceux passant par les reliefs 
ou obstacles que le robot doit contourner. Cependant, une transformee de., 
distance operant par propagation utilisee avec une contrainte statique dans 
le cadre de la navigation terrain pour robot, ne convient pas a la navigation . 
terrain pour aeronef pour lequel la menace presentee par un relief ou un 
obstacle au sol depend du profil vertical de sa trajectolre. 

La presente invention a pour but un procede d'estimation des 
distances des points d'une carte extraite d'une base de donnees d'elevation 
du terrain par rapport a un point de reference employant une transformee de 
distance operant par propagation, avec une contrainte dynamique evolutive 
au cours du temps convenant a la navigation terrain pour un aeronef ayant 
une trajectoire a profil vertical impose. 

Elle a pour objet un procede d'estimation des distances des points 
d'une carte extraite d'une base de donnees d'elevation du terrain, pour un 
mobile soumis a des contraintes dynamiques lui interdisant certaines zones 
de la carte dites zones .interdites.de passage dont la configuration varie en 
fonction du temps de parcours du mobile. La base de donnees d'elevation du 



terrain renferme un ensemble de points reperes par une altitude, une latitude 
et une longitude maillant le terrain devolution du mobile. Le procede met en 
oeuvre une transformee de distance operant par propagation sur I'image 
constitute par !es elements de la base de donnees d'elevation du terrain 
correspondant a la carte, disposes en lignes et colonnes par ordres de 
valeurs de longitude et de latitude. Cette transformee de distance estime les 
distances des differents points de I'image par rapport a un point source place 
a proximite du mobile en appliquant, par balayage, un masque de chanfrein 
aux differents points de I'image. L'estimation de distance d'un point par 
application du masque de chanfrein a ce point dit d'un point but s'effectue en 
repertoriant les differents trajets allant du point but au point source et passant 
par des points du voisinage du point but qui sont recouverts par le masque 
de chanfrein et dont les distances au point source ont ete prealablement 
estimees au cours du meme balayage, en determinant les longueurs des 
drfferents trajets repertories par sommation de la distance affectee au point 
de passage du voisinage et de sa distance au point but extraite du masque 
de chanfrein, en recherchant le trajet le plus court parmi les trajets 
repertories et en adoptant sa longueur comme estimation de la distance du 
po.nt but. Initialement, en debut de balayage, une valeur de distance 
superieure a la plus grande distance mesurable sur I'image est attribute a 
tous les points de I'image sauf pour le point source, origine des mesures de 
distance, auquel est affectee une valeur de distance nulle. Le procede est 
remarquable en ce que les longueurs des trajets repertories, lors de 
rappl.cat.on du masque de chanfrein a un point but, en vue de la recherche 
du trajet le plus court, sont traduites en temps de parcours pour le mobile et 
en ce que les trajets repertories dont les temps de parcours pour le mobile 
sont tels que le point but appartiendrait a une zone interdite de passage au 
moment ou le mobile I'atteindrait, sont exclus de la recherche du trajet le plus 
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Avantageusement, lorsque le mobile est un aeronef ayant un profil 
vertical de vol a respecter determinant revolution de son altitude instantanee 
on associe, aux longueurs des trajets repertoriees, les valeurs previsibles 
des altitudes instantanees qu'aurait I'aeronef en atteignant le point but par 
ces trajets tout en respectant le profil vertical de vol impose, et on elimine de 
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la recherche du trajet le plus court, les trajets repertories associes a des 
valeurs previsibles d'altitude inferieures ou tgales a celle du point but 
donnee par la base de donntes d'tltvation du terrain et augmentee d'une 
marge de protection. 

Avantageusement, lorsque le mobile est un aeronef ayant un profil 
vertical de vol impose, I'estimation de distance optrte par propagation sur 
I'image constitute des elements de la base de donnees d'tltvation du terrain 
correspondant a la carte est doublee d'une estimation de I'altitude previsible 
pour I'aeronef au droit des differents points de I'image en supposant qu'il 
suive le trajet le plus court stlectionne pour I'estimation de distance et qu'il 
respecte le profil vertical de vol impost. 

Avantageusement, lorsque le mobile est un aeronef a profil 
vertical de vol impost et que I'estimation de distance est doublee par une 
estimation de I'altitude previsible de I'aeronef, les altitudes des differents 
points de la carte sont soustraites des estimations des altitudes previsibles 
pour I'aeronef en ces points pour obtenir des ecarts par rapport au sol. 

Avantageusement, lorsque le mobile est un aeronef a profil 
vertical de vol impose et que I'estimation de distance est doublee par une 
estimation de I'altitude previsible de I'aeronef, les altitudes des differents 
points de la carte sont soustraites des estimations des altitudes previsibles 
pour I'aeronef en ces points pour obtenir des ecarts par rapport au sol 
affiches sur la carte en strates de couleurs. 

Avantageusement, la transformte de distance par propagation 
balaie les pixels de I'image constitute des tltments de la base de donntes 
d'tltvation du terrain correspondant a la carte, en plusieurs passes 
successives selon des ordres differents. 

. Avantageusement, la transformte de distance par propagation 
balaie les pixels de I'image constitute des tltments de la base de donntes 
d'tltvation du terrain appartenant k la carte, en plusieurs passes 
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successives selon des ordres differents et de manure repetee jusqu'a ce que 
les estimations de distance obtenues se stabilisent. 

Avantageusement, la transformee de distance par propagation 
balaie les pixels de I'image constitute des elements de la base de donnees 
d'elevation du terrain correspondant a la carte, en plusieurs passes 
successives selon des ordres differents dont I'ordre lexicographique, I'ordre 
lexicographique inverse, I'ordre lexicographique transpose et ' I'ordre 
lexicographique transpose inverse. 

Avantageusement, la transformee de distance par propagation 
balaie les pixels de I'image constitute des elements de la base de donnees 
d'elevation terrain correspondant a la carte, en une serie de quatre passes 
repetee jusqu'a stabilisation des estimations de distances : 

- une premiere passe effectuee ligne par ligne de haut en bas de 
I'image, chaque ligne etant parcourue de gauche a droite, 

- une deuxieme passe effectuee ligne par ligne de bas en haut 
de I'image, chaque ligne ttant parcourue de droite a gauche, 

- une troisieme passe effectuee colonne par colonne de gauche 
a droite de 1'image, chaque colonne etant parcourue de haut 
en bas, et 

- une quatrieme passe effectuee colonne par colonne de droite a 
gauche de I'image, chaque colonne etant parcourue de bas en 
haut. 



Avantageusement, la transformee de distance par propagation 
balaie les pixels de I'image constitute des elements de la base de donnees 
d'elevation du terrain correspondant a la carte, en une serie de huit passes 
repetee jusqu'a stabilisation des estimations de distances. : 

- une premiere passe effectuee ligne par ligne de haut en bas de 
I'image, chaque ligne etant parcourue de gauche a droite, 

- une deuxieme passe effectuee ligne par ligne de bas en haut 
de I'image, chaque ligne etant parcourue de droite a gauche, 
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- une troisieme passe effectuee colonne par colonne de gauche 
a droite de I'image, chaque colonne <§tant parcourue de haut 
en bas, 

- une quatrieme passe effectuee colonne par colonne de droite a 
gauche de I'image, chaque colonne etant parcourue de bas en 
haut, 

- une cinquieme passe effectuee ligne par ligne de haut en bas 
de I'image, chaque ligne etant parcourue de droite a gauche, 

- une sixieme passe effectuee ligne par ligne de bas en haut de 
I'image, chaque ligne etant parcourue de gauche a droite, 

- une septteme passe effectuee colonne par colonne de droite a 
gauche de I'image, chaque colonne etant parcourue de haut en 
bas, et 

- une huitieme passe effectuee colonne par colonne de gauche 

a droite de I'image, chaque colonne etant parcourue de bas en v 
haut. 

D'autres caracteristiques et avantages de I'invention ressortiront. v - 
de la description ci-apres d'un mode de realisation donne a titre d'exemple. i 
Cette description sera faite en regard du dessin dans lequel : 

- une figure 1 represente un exemple de masque de chanfrein, 

- des figures 2a et 2b montrent les cellules du masque de 
chanfrein illustre a la figure 1, qui sont utilisees dans une 
passe de balayage selon I'ordre lexicographique et dans une 
passe de balayage selon I'ordre lexicographique inverse, 

- une figure 3 est un diagramme illustrant les principales etapes 
d'un procede, conforme a I'invention, pour estimer la distance 
d'un point en tenant compte d'une contrainte dynamique au 
cours de Implication d'un masque de chanfrein, 

- une figure 4 est un diagramme illustrant une variante du 
procede d'estimatibn de la distance d'un point montrS a la 
figure 3, et 

. - une figure .5 est un diagramme des principales etapes d'un 
procede, conforme a I'invention, d'estimation, par propagation, 
des distances de I'ensemble des points d'une carte tenant 
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compte d'une contrainte dynamique et mettant en oeuvre un 
precede d'estimation de la distance d'un point tel que ceux 
montres aux figures 3 et 4. 

La distance d'entre deux points d'une surface est la longueur 
minimale de tous les parcours possibles sur la surface partant de Tun des 
points et aboutissant a I'autre. Dans une image formee de pixels repartis 
selon un maillage regulier de lignes, colonnes et diagonales, une 
transformee de distance par propagation estime la distance d'un pixel dit 
pixel "but" par rapport a un pixel dit pixel "source" en construisant 
progressivement, en partant du pixel source, le plus court trajet possible 
suivant le maillage des pixels et aboutissant au pixel but et en s'aidant des 
distances trouvees pour les pixels de I'image deja analyses et d'un tableau 
dit masque de chanfrein repertoriant les valeurs des distances entre un pixel 
et ses proches voisins. 

Comme montre a la figure 1, un masque de chanfrein se presente 
sous la forme d'un tableau avec une disposition de cases reproduisant le 
motif d'un pixel entoure de ses proches voisins. Au centre du motif, une case 
affectee de la valeur 0 repere le pixel pris pour origine des distances 
repertoriees dans le tableau. Autour de cette case centrale, s'agglomerent 
des cases peripheriques remplies de valeurs de distance non nulles et 
reprenant la disposition des pixels du voisinage d'un pixel suppose occuper 
la case centrale. La valeur de distance figurant dans une case peripherique 
est celle de la distance separant un pixel occupant la position de la case 
peripherique concernee, d'un pixel occupant la position de la case centrale 
On remarque que les valeurs de distance se repartissent en cercles 
concentriques. Un premier cercle de quatre cases correspondant aux quatre 
pixels les plus proches du pixel de la case centrale places soit sur la ligne 
soit sur la colonne du pixel de la case centrale sent affectees d'une valeur de 
distance D1 . Un deuxieme cercle de quatre cases correspondant aux quatre 
pixels les plus proches du pixel de la case centrale places en dehors de la 
ligne et de la colonne du pixel de la case centrale sont affectees d'une valeur 
de distance D2. Un troisieme cercle de huit cases correspondant aux huit 
pixels les plus proches du pixel de la case centrale places en dehors de la 
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ligne, de la colonne et des diagonales du . pixel de la case centrale sont 
affectees d'une valeur D3. 

Le masque de chanfrein peut couvrir un voisinage plus ou moins 
etendu du pixel de la case centrale en repertoriant les valeurs des distances 
d'un nornbre plus ou moins important de cercles concentriques de pixels du 
voisinage. II peut etre reduit aux deux premiers cercles formes par les pixels 
du voisinage d'un pixel occupant la case centrale ou etre etendu au-dela des 
trois premiers cercles form§s par les pixels du voisinage du pixel de la case 
centrale mais il est habituel de s'arreter a trois premiers cercles comme celui 
represents a la figure 1. Les valeurs des distances D1, D2, D3 qui 
correspondent a des distances euclidiennes sont exprimees dans une 
echelle autorisant I'emploi de nombres entiers au prix d'une certaine 
approximation. C'est ainsi que G. Borgefors donne a la distance d1 
correspondant a un echelon en abscisse x ou en ordonnee y la valeur 5, a la, 
distance d2 correspondant a la racine de la somme des carres des echelons 
en abscisse et ordonnee ^x 2 + y 2 la valeur 7 qui est une approximation de * 
5V2 , et a la distance d3 la valeur 1 1 qui est une approximation de 5<Js . 

La construction progressive du plus court trajet possible allant a 
un pixel but, en partant d'un pixel source et en suivant le maillage des pixels > 
se fait par un balayage regulier des pixels de Pimage au moyen du masque: 
de chanfrein. Initialement, les pixels de Pimage se voient affecter une valeur 
de distance infinie, en fait un nombre suffisamment eleve pour depasser. 
toutes les valeurs des distances mesurables dans Pimage, a Pexception du 
pixel source qui se voit affecter une valeur de distance nulle. Puis les valeurs 
initiates de distance affectees aux points but sont mises a jour au cours du 
balayage de Pimage par le masque de chanfrein, une mise a jour consistant 
a remplacer une valeur de distance attribute a un point but, par une nouvelle 
valeur moindre resultant d'une estimation de distance faite a Poccasion d'une 
nouvelle application du masque de chanfrein au point but consider^. 

Une estimation de distance par application du masque de 
chanfrein a un pixel but consiste a repertorier tous les trajets allant de ce 
pixel but au pixel source et passant par un pixel du voisinage du pixel but 
dont la distance a deja ete estimee au cours du meme balayage, a 
. . rechercher parmi les trajets repertories, le ou les trajets les plus courts et a 
adopter la longueur du ou des trajets les plus courts comme estimation de 
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distance. Cela se fait en placant le pixel but dont on veut estimer la distance 
dans la case centrale du masque de chanfrein, en selectionnant les cases 
per.pher.ques du masque de chanfrein correspondant a des pixels du 
vo.s.nage dont la distance vient d'etre mise a jour, en calculant les longueurs 
des trajets les plus courts reliant le pixel a mettre a jour au pixel source en 
passant par un des pixels selectionnes du voisinage, par addition de la 
valeur de distance affectee au pixel du voisinage concerne et de la valeur de 
distance donnee par le masque de chanfrein, et a adopter, comme 
est.mat,on de distance, le minimum des valeurs de longueur de trajet 
obtenues et de I'ancienne valeur de distance affectee au pixel en cours 
d'analyse. 

L'ordre de balayage des pixels de I'image influe sur la fiabilite des 
estimations de distance et de leurs mises a jour car les trajets pris en compte 
en dependent. En fait, il est soumis a une contrainte de regularite qui fait que 
s. les p,xels de I'image sont reperes selon l'ordre lexicographique (pixels 
classes dans un ordre croissant ligne par ligne en partant du haut de I'image 
et en progressant vers le bas de I'image, et de gauche a droite au sein d'une 
l«gne), et si un pixel p a ete analyse avant un pixel q alors un pixel p +x doit 
etre analyses avant le pixel q+ x. Les ordres lexicographique (balayage des 
p.xels de I'image ligne par ligne de haut en bas et, au sein d'une ligne de 
gauche a droite), lexicographique inverse (balayage des pixels de I'image 
.gne par ligne de bas en haut et, au sein d'une ligne, de droite a gauche) 
lex.cographique transpose (balayage des pixels de I'image colonne par 
colonne de gauche a droite et, au sein d'une colonne, de haut en bas) 
lexicographique transpose inverse (balayage des pixels par colonnes de 
dro.te a gauche et au sein d'une colonne de bas en haut) satisfont cette 
cond,t,on de regularite et plus generalement tous les balayages dans 
lesquels les lignes et colonnes, ou les diagonals sont balayees de droite a 
gauche ou de gauche a droite. G. Borgefors preconise un double balayage 
des p,xels de I'image, une fois dans l'ordre lexicographique et une autre dans 
I ordre lexicographique inverse. 

La figure 2a montre, dans le cas d'une passe de balayage selon 
I ordre lex.cographique allant du coin superieur gauche au coin inferieur droit 
de I, mage, les cases du masque de chanfrein de la figure 1 utilisees pour 
repertoner les trajets allant d'un pixel but place sur la case centrale (case 
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indexee par 0) au pixel source en passant par un pixel du voisinage dont la 
distance a deja fait I'objet d'une estimation au cours du meme balayage. Ces 
cases sont au nombre de huit, disposees dans la partie superieure gauche 
du masque de chanfrein. II y a done huit trajets repertories pour la recherche 
du plus court dont la longueur est prise pour estimation de la distance. 

La figure 2b montre, dans le cas d'une passe de balayage selon 
I'ordre lexicographique inverse allant du coin inferieur droit au coin superieur 
gauche de I'image, les cases du masque de chanfrein de la figure 1 utilisees 
pour repertorier les. trajets allant d'un pixel but place sur la case centrale 
(case indexee par 0) au pixel source en passant par un pixel du voisinage 
dont la distance a deja fait I'objet d'une estimation au cours du meme 
balayage. Ces cases sont complementaires de celles de la figure 2a. Elles 
sont egalement au nombre de huit mais disposees dans la partie inferieure 
droite du masque de chanfrein. II y a done encore huit trajets repertories pour 
la recherche du plus court dont la longueur est prise pour estimation de la, 
distance. 

La transformee de distance par propagation dont le principe vient... 
d'etre rappele sommairement a ate congue a I'origine pour I'analyse du v 
positionnement d'objets dans une image mais elle n'a pas tarde a etre- 
appliquee a I'estimation des distances sur une carte du relief extraite d'une i 
base de donnees d'elevation du terrain a maillage regulier de la surface 
terrestre. En effet, une telle carte ne dispose pas explicitement d'une 
m&rique puisqu'elle est tracee a partir des altitudes des points du maillage 
de la base de donnees d'elevation du terrain de la zone representee. Dans 
ce cadre, la transformee de distance par propagation est appliquee a une 
image dont les pixels sont les elements de la base de donnees d'elevation du 
terrain appartenant a la carte, e'est-a-dire, des valeurs d'altitude associees 
aux coordonnees geographiques latitude, longitude des nceuds du maillage 
ou elles ont ete mesurees, classes, comme sur la carte, par latitude et par 
longitude croissantes ou d6croissantes selon un tableau a deux dimensions 
de coordonnees latitude et longitude. 

Pour une navigation terrain de mobiles tels que des robots, la 
transformee de distance par propagation, est utilisee pour estimer les 
distances des points de la carte d'un terrain d'evolution extraite d'une base 
de donnees d'elevation du terrain par. rapport a la position du mobile ou une 
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position proche. Dans ce cas, il est connu de tenir compte de contraintes 
statiques constitutes par des zones de la carte infranchissables par le 
mobile en raison de leurs configurations accidentees. Pour ce faire, un 
marqueur de zone interdite est associe aux elements de la base de donnees 
d'elevation du terrain figurant dans la carte. II signale, lorsqu'il est active, une 
zone infranchissable ou interdite et inhibe toute mise a jour autre qu'une 
initialisation, de I'estimation de distance faite par la transformee de distance 
par propagation pour Pelement pixel considere. 

Dans le cas d'un aeronef, les zones infranchissables evoluent en 
fonction du profll vertical impose a sa trajectoire si bien qu'une estimation de 
distance sous contraintes statiques au moyen d'une transformee de distance 
par propagation n'est pas satisfaisante. 

On propose de prendre en compte, dans la definition des zones 
>nterd.tes de passage, Paltitude previsible de I'aeronef a chaque point but 
dont la distance est en cours d'estimation. Cette altitude previsible qui 
depend bien evidemment du trajet emprunte, est celle de I'aeronef apres 
suivi du trajet adopte pour la mesure de distance, (.'estimation de cette 
altrtude previsible de I'aeronef en un point but, se fait par propagation au 
cours du balayage de f'image par le masque de chanfrein d'une maniere 
analogue a I'estimation de distance. Pour chaque trajet repertorie allant d'un 
point but au point source en passant par un point du voisinage du point but 
dont la distance au point source et ('altitude previsible de I'aeronef ont deja 
ete estimees au cours du meme balayage, ('altitude previsible de I'aeronef 
est deduite de la longueur du trajet et du profil vertical impose a la trajectoire 
de I'aeronef. Cette altitude previsible, estimee pour chaque trajet repertorie 
allant d'un point but dont la distance est en cours d'estimation a un point 
source place a proximite de la position de I'aeronef, est utilisee comme un 
cntere de selection des trajets pris en compte dans I'estimation distance Si 
elle est inferieure ou egale a Paltitude du point but figurant dans la base de 
donnees d'elevation du terrain majoree d'une marge de securite, le trajet 
repertorie auquel elle est associee est ecarte et ne participe pas a la 
selection du plus court trajet. Une fois la selection du plus court trajet 
effectuee, sa longueur est prise pour distance du point but et Paltitude 
previsible pour I'aeronef qui lui est associee est egalement retenue pour 
Paltitude de I'aeronef au point but. 
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La figure 3 illustre les principles etapes du. traitement.effectue 
lors de ^application du masque de chanfrein a un point but P fJ pour estimer 
sa distance pour un aeronef ayant un profil vertical de trajectoire impose. Le 
point but considere P u est place dans la case centrale du masque de 
chanfrein. Pour chaque point voisin P v qui entre dans les cases du masque 
de chanfrein et dont la distance a deja ete estimee au cours du meme 
balayage, le traitement consiste a : 

- lire la distance estimee D v du point voisin Pv (etape 30), 

- lire I'altitude Ag du point but Pij dans la base de donnees 
d'elevation du terrain (etape 31), 

- lire le coefficient Cxy du masque de chanfrein correspondant a 
la case occupee par le point voisin P v (etape 32), 

- calculer la distance propagee D P correspondant a la somme de 
la distance estimee D v du point voisin Pv et du coefficient Cxy 
affecte a la case du masque de chanfrein occupee par le point, 
voisin Pv : 

D^Dy+Cxy 

(etape 33), 

- calculer I'altitude previsible A P de I'aeronef apres . 
franchissement de la distance D P directement a partir de la -J 
distance D P si le profil vertical impose a la trajectoire de 
I'aeronef est defini en fonction de la distance parcourue PV(D P ) 
et prend implicitement en compte le temps de parcours ou 
indirectement par I'intermediaire du temps de parcours si le 
profil vertical impose a la trajectoire de I'aeronef est defini par 
une Vitesse de changement d'altitude (etape 34), 

- comparer I'altitude previsible A P obtenue avec celle Ay du point 
but Pij tiree de la base de donnees d'elevation du terrain 
augmentee d'une marge de securite A (elape 35), 

- eliminer. la distance propagee D P si I'altitude previsible A P est 
inferieure ou egale a celle Ay du point but Py tiree de la base 
de donnees d'elevation du terrain et majoree par la marge de 

. securite. A (etape 36), 

- si I'altitude previsible A P est supeiieure a celle Ay du point but 
Py majoree par la marge de securite A, lire la distance Dg deja 
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affectee au point but consider^ Py ( etape 37) et la comparer a 
la distance propagee D P (etape 38), 

- eliminer la distance propagee D P si elle est superieure ou 
egale a la distance Dy deja affectee au point but consider^ p, Jf 
et 

- remplacer la distance Dy deja affectee au point but consider^ 
Pi.j, par la distance propagee D P si cette derniere est inferieure 
(etape 39). 

La figure 4 illustre les principals etapes d'une variante du 
traitement effectue lors de ('application du masque de chanfrein a un point 
but P u pour estimer sa distance pour un aeronef ayant un profil vertical de 
trajectoire impose. 

Cette variante differe dans la maniere d'elaborer I'altitude 
previsible A P de I'aeronef. Elle suppose que I'altitude previsible pour I'aeronef 
en chaque point de la base de donnees d'elevation du terrain calculee en 
fonction du profil vertical impose a sa trajectoire et a partir de la longueur du 
trajet selectionne pour la mesure de distance n'est pas consideree comme 
une variable passagere, ce que permet le traitement decrit relativement a la 
figure 3,. mais est memorisee, au meme titre que I'estimation de distance. 
Dans cette variante, les etapes 30, 31 de lectures de la distance estimee D v 
du point voisin P v et de I'altitude A,j du point but Py dans la base de donnees 
d'elevation du terrain sont completes par une etape 40 de lecture de 
I'altitude previsible A PV pour I'aeronef au point voisin P v> et le calcul de 
I'altitude previsible A P se fait (etape 34') par sommation de I'altitude previsible 
A PV au point voisin P v et de la variation d'altitude sur la distance separant le 
point voisin P v du point but P g due au profil vertical impose a la trajectoire de 
I'aeronef. 

La memorisation des altitudes previsibles pour i'aeronef lorsqu'il 
atteint les differents points de la carte qui lui sont accessibles permet 
d'etablir, en leur soustrayant les altitudes des points de la carte tirees de la 
base de donnees d'elevation du terrain, une carte des possibilites maximales 
de hauteurs de survol de I'aeronef representant les ecarts previsibles par 
rapport au terrain en strates de couleurs. Une telle carte facilite pour 
I'equipage de I'aeronef le choix d'une trajectoire realiste presentant la 
meilleure garde au sol. 
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Comme indique precedemment, I'estimation des distances des 
differents points de la carte se fait en appliquant un traitement par masque 
de chanfrein tel que ceux qui viennent d'etre decrits relativement aux figures 
3 et 4, a I'ensemble des pixels de I'image formee par les elements de la base 
de donnees d'elevation du terrain appartenant a la carte, pris 
successivement selon un balayage regulier comportant un minimum de deux 
passes realisees dans des ordres inverses. 

La figure 5 illustre les principales etapes d'un exemple de 
processus global permettant I'estimation des distances de I'ensemble des 
points d'une carte de relief pour un mobile soumis a des contraintes 
dynamiques. 

La premiere etape 50 du processus est une initialisation des 
distances affectees aux differents points de la carte consideres comme les 
pixels d'une image. Cette initialisation des distances consiste, comme. 
indique precedemment, a attribuer une valeur de distance infinie, a tout le . 
moins superieure a la plus grande distance mesurable sur la carte, pour tous ; 
les points de la carte consideres comme des points but, a I'exception d'un,, 
seul considere comme la source de toutes les distances auquel est attribu6 
une valeur de distance nulle. Ce point source est choisi a proximity de la . 
position instantanee du mobile sur la carte. 

Les etapes suivantes 51 a 54 sont des passes d'un balayage .; 
regulier au cours desquelles le masque de chanfrein est appliqu6 
successivement et a plusieurs reprises a tous les points de la carte 
consideres comme les pixels d'une image, I'application du masque de 
chanfrein a un point de la carte donnant a une estimation de la distance de 
ce point par rapport au point source, par execution d'un des traitements 
decrits relativement a la figure 3 ou a la figure 4. 

La premiere passe de balayage (etape 51) se fait dans I'ordre 
lexicographique, les pixels de I'image etant analyses ligne par ligne du haut 
vers le bas de I'image et de gauche a droite au sein d'une m§me ligne. La 
deuxieme passe de balayage (etape 52) se fait dans I'ordre lexicographique 
inverse, les pixels de I'image etant toujours analyses ligne par ligne mais du 
bas vers le haut. de I'image et de droite a gauche au sein d'une ligne. La 
troisieme passe de balayage (6tape 53) se fait dans I'ordre lexicographique 
transpose, les pixels de I'image etant analyses colonne par colonne de la 



16 



gauche vers la droite de I'image et de haut en bas au sein d'une meme 
colonne. La quatrieme passe de balayage (<§tape 54) se fait dans I'ordre 
lexicographique transpose inverse, les pixels de I'image etant analyses 
colonne par colonne mais de la droite vers la gauche de I'image et de bas en 
haut au sein d'une meme colonne. 

Ces quatre passes ( etapes 51 a 54) sont repetees tant que 
I'image de distance obtenue change. Pour ce faire, le contenu de I'image de 
distance obtenu est memorise (<§tape 56) apres chaque serie de quatre 
passes ( etapes 51 a 54) et compare avec le contenu de I'image de distance 
obtenu a la s§rie precedente (tope 55), la boucle n'etant brisee que lorsque 
la comparaison montre que le contenu de I'image distance ne varie plus. 

En theorie, deux passes de balayage selon I'ordre lexicographique 
et I'ordre lexicographique inverse peuvent suffire. Cependant la presence de 
zones interdites de passage de forme concave peut provoquer, dans le 
phenomene de propagation des distances, des angles morts renfermant des 
pixels pour lequel ('application du masque de chanfrein ne donne pas 
d'estimation distance. Pour diminuer ce risque d'angle mort, il y a lieu de faire 
vaner la direction du phenomene de propagation distance en faisant varier la 
d.rect.on du balayage d'ou le doublement des passes avec une transposition 
des ordres de balayage correspondant a une rotation de I'image de 90° Pour 
une encore meilleure elimination des angles morts, on peut proceder a des 
series de huit passes : 

- une premiere passe effectuee ligne par ligne de haut en bas de 
I'image, chaque ligne etant parcourue de gauche a droite, 

- une deuxieme passe effectuee ligne par ligne de bas en haut 
de I'image, chaque ligne etant parcourue de droite a gauche, 

- une troisieme passe effectuee colonne par colonne de gauche 
a droite de I'image, chaque colonne etant parcourue de haut 
en bas, 

- une quatrieme passe effectuee colonne par colonne de droite a 
gauche de I'image, chaque colonne etant parcourue de bas en 
haut, 

- une cinquieme passe effectuee ligne par ligne de haut en bas 
de 1'image, chaque ligne etant parcourue de droite a gauche, 
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- une sixieme passe effectuee ligne par ligne de bas en .haut de 
I'image, chaque ligne etant parcourue de gauche a droite, 

- une septieme passe effectuee colonne par colonne de droite a 
gauche de I'image, chaque colonne etant parcourue de haut en 
bas, et 

- une huitieme passe effectuee colonne par colonne de gauche 
a droite de I'image, chaque colonne etant parcourue de bas en 
haut. 

II est possible d'introduire dans la serie de passes de balayage 
d'autres types de passes de balayage deduits des passes precedentes en 
faisant jouer par les diagonales de I'image, les roles joues precedemment 
par les lignes et colonnes de I'image. Cela revient a appliquer les passes de 
balayage decrites precedemment a une image tournee de 45°. De maniere 
generate, plus les passes d'une serie sont variees plus le risque d'angle mort 
diminue. 
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REVENDICATIONS 



1. Precede d'estimation des distances des points d'une carte 
extraite d'une base de donnees d'elevation du terrain, pour un mobile soumis 
a des contraintes dynamiques lui interdisant certaines zones de la carte dites 
zones interdites de passage dont la configuration varie en fonction du temps 
de parcours du mobile ; la base de donnees d'elevation du terrain renfermant 
un ensemble de points reperes par une altitude, une latitude et une longitude 
maillant le terrain devolution du mobile ; ledit precede mettant en ceuvre une 
transformee de distance operant par propagation sur I'image constitute par 
les elements de la base de donnees d'elevation du terrain correspondant a la 
carte, disposes en lignes et colonnes par ordres de valeurs de longitude et 
de latitude ; la transformee de distance estimant les distances des differents 
points de I'image par rapport a un point source place a proximite du mobile, 
en appliquant, par balayage, un masque de chanfrein aux differents points de 
I'image ; I'estimation de distance d'un point, par application du masque de 
chanfrein a ce point dit point but s'effectuant en repertoriant les differents 
trajets allant du point but au point source et passant par des points du 
voisinage du point but qui sont couverts par le masque de chanfrein et dont 
les distances au point source ont ete prealablement estimees au cours du 
meme balayage, en determinant les longueurs des differents trajets 
repertories par sommation de la distance affectee au point de passage du 
voisinage et de sa distance au point but extraite du masque de chanfrein en 
recherchant le trajet le plus court parmi les trajets repertories et en adoptant 
sa longueur comme estimation de la distance du point but ; une valeur de 
distance superieure a la plus grande distance mesurable sur I'image etant 
m.tialement attribute, en debut de balayage, a tous les points de I'image sauf 
au point source, origine des mesures de distance, auquel est affectee une 
valeur de distance nulle ; ledit precede etant caracterise en ce que les 
longueurs des trajets repertories, lors de ^application du masque de 
chanfrein a un point but, en vue de la recherche du trajet le plus court sont 
tradu.tes en temps de parcours pour le mobile et en ce que les trajets 
repertories dont les temps de parcours pour le mobile sont tels que le point 
but appartiendrait a une zone interdite de passage au moment ou le mobile 
I'atteindrait, sont exclus de la recherche du trajet le plus court. 
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2. Precede selon la revendication 1 , applique a un aeronef ayant 
un profil vertical de vol a respecter determinant revolution de son altitude 
instantanee, caracterise en ce que, aux longueurs des trajets repertoriees 
lors de ^'application du masque de chanfrein a un point but, sont associees 
les valeurs previsibles des altitudes instantanees qu'aurait I'aeronef en 
atteignant le point but par ces trajets tout en respectant le profil vertical de 
vol impose, et en ce que les trajets repertories associes a des valeurs 
previsibles d'altitude inferieures ou egales a celle du point but donnee par la 
base de donnees d'elevation du terrain et augmentee d'une marge de 
protection sont exclus de la recherche du trajet le plus court,. 

3. Procede selon la revendication 2, caracterise en ce que 
l'estimation de distance operee par propagation sur I'image constitute des 
elements de la base de donnees d'elevation du terrain correspondant a la 
carte est doublee d'une estimation de I'altitude previsible pour I'aeronef au 
droit des differents points de I'image en supposant qu'il suive le trajet le plus 
court selectionne pour l'estimation de distance et qu'il respecte le profil 
vertical de vol impose. 

4. Procede selon la revendication 3, caracterise en ce que les 
altitudes des differents points de la carte sont soustraites des estimations 
des altitudes previsibles pour I'aeronef en ces points pour obtenir des ecarts 
par rapport au sol. 

5. Procede selon la revendication 4, caracterise en ce que les 
ecarts par rapport au sol sont affiches sur la carte en strates de couleurs. 

6. Procede selon la revendication 1 , caracterise en ce que la 
transformee de distance par propagation balaie les pixels de I'image 
constitute des elements de la base de donnees d'elevation du terrain 
correspondant a la carte, en plusieurs passes successives selon des ordres 
differents. 
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7. Procedt selon la revendication 6, caracterise en ce que la 
transformee de distance par propagation balaie les pixels de I'image 
constitute des elements de la base de donnees d'elevation du terrain 
appartenant a la carte, en plusieurs passes successives selon des ordres 
differents et de maniere repetee jusqu'a ce que les estimations de distance 
obtenues se stabilisent. 

8. Procede selon la revendication 6, caracterise en ce que la 
transformee de distance par propagation balaie les pixels de I'image 
constitute des elements de la base de donnees d'elevation du terrain 
correspondant a la carte, en plusieurs passes successives selon des ordres 
differents dont I'ordre lexicographique, I'ordre lexicographique inverse, I'ordre 
lexicographique transpose et I'ordre lexicographique transpose inverse. 

9. Procede selon la revendication 6, caracterise en ce que la 
transformee de distance par propagation balaie les pixels de I'image 
constitute des elements de la base de donnees d'elevation du terrain 
correspondant a la carte, en une serie de quatre passes repetee jusqu'a 
stabilisation des estimations de distances. : 

- une premiere passe effectuee ligne par ligne de haut en bas de 
I'image, chaque ligne etant parcourue de gauche a droite, 

- une deuxieme passe effectuee ligne par ligne de bas en haut 
de I'image, chaque ligne etant parcourue de droite a gauche, 

- une troisieme passe effectuee colonne par colonne de gauche 
a droite de I'image, chaque colonne etant parcourue de haut 
en bas, et 

- une quatrieme passe effectuee colonne par colonne de droite a 
gauche de I'image, chaque colonne etant parcourue de bas en 
haut. 

10. Procede selon la revendication 6, caracterise en ce que la 
transformee de distance par propagation balaie les pixels de I'image 
constituee des elements de la base de donnees d'elevation du terrain 
correspondant a la carte, en une serie de huit passes repetee jusqu'a 
stabilisation des estimations de distances. : 
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- une premiere passe effectuee ligne par ligne de haut en bas de 
I'image, chaque ligne etant parcourue de gauche a droite, 

- une deuxieme passe effectuee ligne par ligne de bas en haut 
de I'image, chaque ligne etant parcourue de droite a gauche, 

- une troisieme passe effectuee colonne par colonne de gauche 
a droite de I'image, chaque colonne etant parcourue de haut 
en bas, 

- une quatrieme passe effectuee colonne par colonne de droite a 
gauche de I'image, chaque colonne etant parcourue de bas en 
haut, 

- une cinquieme passe effectuee ligne par ligne de haut en bas 
de I'image, chaque ligne etant parcourue de droite a gauche, 

- une sixieme passe effectuee ligne par ligne de bas en haut de 
I'image, chaque ligne etant parcourue de gauche a droite, 

- une septieme passe effectuee colonne par colonne de droite a 
gauche de I'image, chaque colonne etant parcourue de haut en 
bas, et 

- une huitieme passe effectuee colonne par colonne de gauche 
a droite de I'image, chaque colonne etant parcourue de bas en 
haut. 

i. 

1 1 . Precede selon la revendication 6, caracterise en ce que la 
transformee de distance par propagation balaie les pixels de I'image 
constitute des elements de la base de donnees d'elevation du terrain 
appartenant a la carte, en plusieurs passes successives selon des ordres 
differents dont certains consistent en un balayage de I'image par diagonales, 
d'un bord a I'autre et, au sein d'une diagonale, d'une extremite a I'autre. 
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